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We  consider  the  motion  of  vortex  blotches,  a two  dimensional 
incompressible  flow  with  zero  vorticity  except  inside  certain 
closed  curves.  The  streamf unction  of  such  a flow  is1 


<p  (x,t) 


U)  (x;  t)  ln|  X-Xj^  | 


(1) 


where  t for  time  and  x denotes  a two-dimensional  position  vector, 
and  the  integral  is  over  areas  where  the  vorticity  u(x,t)  is  different 
from  zero.  We  consider,  for  simplicity,  such  areas  as  a single 
closed  curve  C of  constant  vorticity  within.  The  velocity  field 
is  then  given  by 


u(x,t)  = ” (j?  Mxj.)  ln|x-x1|dJl1  (2) 

A 

where  t is  a unit  tangent  counter-clockwise  along  the  curve.  Since 
fluid  particles  on  the  curve  C remain  there  all  the  time,  the  motion 
of  a vortex  blotch  is  completely  described  by  the  motion  of  these 
particles.  Denoting  the  position  vector  of  points  on  curve  C 
by  x * [x (t,s) ,y (t,s) ] , where  s is  a parametric  coordinate  on  C, 
we  have  the  equation  of  motion  as 


n = “ Ir £ ^(?i)lnir?iidsi 

where  x1  = [x  (t,  s^  ,y  (t,  s^  ] . 


(3) 


2. 


It  can  be  shown,  for  any  C that  the  following  integrals 
along  C are  constant  in  time  : 


Area  within  = | (x  dy  ••  y dx), 
Centroid  = (■|^J’x^dy,  - -^(^J"y2  dx). 


Moment  of  inertia 


ia  = i (j>r  (x3dy  - y3dx)  . 


In  this  note  we  concern  ourselves  with  the  question  of  in- 
variant curves,  i.e.  those  curves  which  do  not  change  their 
shape  in  time.  The  motion  of  an  invariant  curve  must  be  a com- 
bination of  solid  translation,  solid  rotation,  and  sliding  along 
the  curve.  The  translational  motion  is  prohibited  by  the  in- 
variant of  centroid  as  listed  herein.  The  condition  for  invariant 
curves  can  be  expressed  as  in  a certain  rotating  frame  of  refer- 
ence, the  normal  velocity  of  every  point  or  the  curve  should 
vanish.  Mathematically  this  becomes,  after  using  (3)  in  a 
rotating  frame  of  reference  about  the  centroid  of  C with  an 
angular  speed  ft  , 

ft  t (x ) • x = ~ (jjr  In  Jx-x^  ri(x)  •t(x1)ds1  (4) 

when  n is  an  outward  unit  normal. 

In  terms  of  the  parametric  representation  of  the  curve 


x = x(s)  , y = y (s)  and  xx  = x(sx)  , y1  = y(s1) 


we  have 


X |^(x2+y2)  = ln|  Jj-Xjl  (yxj-xy^ds.,^  (5) 

where  a dot  denotes  differentiation  with  its  argument  and  *=  ftir/t 


3. 


. 


Equation  (5)  takes  the  form  of  the  eigenvalue  problems;  all  the 
eigenfunctions  represent  the  possible  invariant  curves,  and  the 
corresponding  eigenvalue  give  the  speed  of  the  rotating  frame 
wherein  the  curve  remains  fixed  in  time. 

Equation  (5)  is  invariant  under  a uniform  stretching,  and 
a solid  rotation.  Therefore  the  size  and  the  orientation  of 
the  figure  are  immaterial  as  one  would  expect.  It  is, however, 
not  a translational  invariant.  The  origin  (x  = 0,  y = 0)  has 
been  taken  at  the  centroid  of  the  figure.  This  introduces  two 
integral  constraints 


* 


and 


y dx  = 0 


(6) 


x dy  = 0 
to  the  solution  of  (5) . 

We  are  looking  for  periodic  solutions  of  x(s)  and  y(s).  It 
is  natural  to  express  these  in  Fourier  series. 


x(6)  = £ (a  cos  n0  + b sin  n0)  , 

1X1  XI 


y(0)  = 7 (a  cos  n©  + 3^  sin  n0)  , 

**  i n n 

n=l 


(7) 


where  we  choose  0 instead  of  s,  ranges  from  0 to  2F  as  the 
parameter.  The  logarithmic  term  under  the  integral  sign  of 
(5)  make  the  integral  quite  difficult  to  evaluate.  However, 
if  one  takes  only  one  value  of  n in  (7) , one  can  solve  (5) 
explicitly.  Now  for  a single  value  of  n,  Eq.  (7)  represents 
an  ellipse.  Realizing  that  (5)  is  Indifferent  to  the  orientation 
we  take  (7)  for  a single  value  of  n to  be  simply 
x = a cos  0 , y = b sin  0 . 


This  solves  (5)  provided  A = irab/(a+b)'  • The  speed  of  rotating 
frame  of  reference  is  then  given  by 


fl 


= wab 


2 ' 


(B) 


(a+b) 

which  checks  with  Kirchhoff's  solution  of  elliptical  vortex1. 

In  a polar  coordinate,  a curve  is  represented  by  r = r (6) , 
and  Eq.  (5)  can  be  put  into  the 

2ir 


dr  _ 1 dr  f ..  2,.. 

2Ar  d©  2r  d«  J0  * 40 

r . _1  _1  f2lT  d(j)  rn+l  } cos  (n_1}  (9) 

de  n=2  n2-l  r11-1  ' 0 


where  we  have  used  a series  expansion  for  the  logarithmic  term  in 
(5)  as  follows: 

2 . J2  „ ..,1/2  _ , V /Ii^n  cos_nM  2 (10 


y z _ , j./ z - r j 

ln[r  + rf  - 2rr.cos  y]  = In  r - \ (—)  — r » 

1 1 n=l  L 


with  r = r (0)  , r^  = r(4>)  and  y = 0-<J> . 


5 


We  next  consider  a perturbation  solution  to  (9)  corresponding 

to  some  small  distortion  of  a circular  vortex.  It  is  found  that 

the  solution  can  be  represented  as  follows: 

2 3 4 

r(Q)  = a0+eai  cos  m9+e  a2  cos  2m8+e  a^  cos  3m0+O(e  ) (11) 

A = AQ+e2A2  (12) 

where  m is  an  integer  taking  to  be  greater  than  one  in  order  for 
(11)  to  satisfy  (6) . We  consider  that  e <<  1 so  that  (11) 
represents  some  small  distortion  of  a circular  vortex. 
Substituting  (11)  and  (12)  into  (9) , and  equating  the  coefficient 

3 

of  each  power  of  e (up  to  e in  the  present  case)  to  zero.  We 
obtain 


Xo  - 

7t  (m-1) 

2m 

2 

(13) 

X2  = 

ir (m-1)  ,al% 
~ 4 (a 

(14) 

o 

, 

a 

2 

2m- 1 ,al. 

4 ' 

(15) 

o 

o 

3 

!i  - 

( 3m- 1 ) (m-1) 

al 

<iT> 

(16) 

a 

8 

o 

o 

The  first 

expression 

above 

checks  with  the  calculation  for 

infinitesimal  perturbation 

of  a circular  vortex  as  given  in 

Lamb1.  The  approximate  invariant  as  given  by  (11)  and  (12)  can 
be  improved  by  adding  higher  terms  in  e.  A systematic  procedure 
for  higher  order  invariant  will  be  presented  elsewhere. 
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